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We study the frequency dependence of the effective electromagnetic parameters of left-handed and related 
meta-materials of the split ring resonator and wire type. We show that the reduced translational symmetry 
(periodic structure) inherent to these meta-materials influences their effective electromagnetic response. To 
anticipate this periodicity, we formulate a periodic effective medium model which enables us to distinguish 
the resonant behavior of electromagnetic parameters from effects of the periodicity of the structure. We use 
this model for the analysis of numerical data for the transmission and reflection of periodic arrays of split ring 
resonators, thin metallic wires, cut wires as well as the left-handed structures. The present method enables us 
to identify the origin of the previously observed resonance/anti-resonance coupling as well as the occurrence of 
negative imaginary parts in the effective permittivities and permeabilities of those materials. Our analysis shows 
that the periodicity of the structure can be neglected only for the wavelength of the electromagnetic wave larger 
than 30 space periods of the investigated structure. 

PACS numbers: 41.20.Jb, 42.25.Bs, 42.70.Qs, 73.20.Mf 



I. INTRODUCTION 



permeability shows a resonant behavior. 



Recent progress in studies of left-handed meta-materials 
(LHM)' confirmed that the fabrication of structures with neg- 
ative effective permittivity and permeability, and their appli- 
cation in technical praxis is possible. The most promising 
structures are based on the combination of periodic arrays of 
metallic split ring resonators (SRR) and thin metallic wires, a 
design proposed theoretically by Pendry et alwi^ and experi- 
mentally verified by Smith et al^^. 

It is assumed that in a well defined frequency interval both 
effective permittivity and permeability of LHM are simulta- 
neously negative. Consequently, also the refractive index is 
negative^. This theoretical prediction was supported experi- 
mentally by measurements of the transmission of the electro- 
magnetic (EM) wave through the LHM: A transmission peak 
was observed in the frequency region where the LH band is 
expected^-^. Negativeness of the index of refraction was ver- 
ified experimentally by the Snell's law experiment^ and con- 
firmed later by other experiments^ '^. Numerical simulations 
were performed which also observed a transmission peak in 
the resonant frequency interval" '^ '^. Effective electromag- 
netic parameters were calculatedii by comparison of numer- 
ically obtained transmission and reflection amplitudes of the 
LHM with theoretical formulas for a homogeneous slab. The 
obtained results confirmed that the refractive index of the 
LHM is indeed negative in the resonant frequency interval. 
Moreover, the obtained frequency dependence of the effective 
permittivity and permeability also agreed qualitatively with 
theoretical predictions. In particular, the effective magnetic 
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typical for lattice of SRR"* in the vicinity of the magnetic reso- 
nance frequency LOm- The effective permittivity is determined 
by the electric response of the array of thin wires. 
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and is negative if the frequency is smaller than the plasma 
frequency ujp. Transmission data, obtained using either the 
transfer matrix method'^ or commercial software'^ was ana- 
lyzed to find the dependence of the resonance frequency uJm 
on the structural parameters of the SRR and on the design of 
the unit cell of the LHM structure. 

Further progress in numerical methods brought more ac- 
curate data and strong evidence that the effective parameters 
of the LHM differ considerably from the theoretical predic- 
tion (^12}. Although the main properties - resonant behavior 
of the magnetic permeability at cj,,, and negativeness of the 
effective permittivity - are clearly visible in the data, the ef- 
fective medium picture is spoiled by partially very significant 
anomalies: 

Resonance/anti-resonance coupling. We expect the electric 
and magnetic response of the discussed meta-materials to be 
independent from each other However, whenever there is 
a resonance in Re we simultaneously observe an anti- 
resonant behavior in Re giiii2i22i^iiS and vice versa2^. 

The anti-resonant structures in the real part are accompanied 
by an negative imaginary part^'i^^. 
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Misshapen, truncated resonances. The divergence in Re n ap- 
pears to be cut-off at the edges of the first Brillouin zone and, 
in particular, the negative regions of the magnetic resonance 
in /i and cut-wire resonance in e do not return from large nega- 
tive real part but seem to saturate in a rather shallow behavior 
The corresponding absorption peak in the imaginary parts is 
misshapen and highly asymmetric tooS. 

Discrepancy between n and z about the positions of the res- 
onances. We expect the peaks (or zeros) in the index of re- 
fraction and the impedance to appear exactly at the resonance 
frequency. From the simulations, however, we find different 
frequencies from n and z, respectively. This leads, for in- 
stance, to an "internal structure" of the magnetic resonance 
as shown in Figs. |5] and |6] This structure can not easily be 
explained within the assumed effective medium pictura^i 

Additional spectral structures. Apart from structures around 
the anticipated contributions of the meta-material's con- 
stituents, we observe a lot of additional structure, especially 
at higher frequency, which can not be accounted for 

The above described observation, especially the negative- 
ness of the imaginary part of effective permittivity or perme- 
ability, raised objections2i2^ of other groups. Efros^ argued 
that the LHM can not be approximated by a homogeneous 
system because of the periodicity of the meta-material'^-^*. 

In this paper, we show that the observed artifacts in the ho- 
mogeneous effective approximation are quite generic. They 
are given by the periodic structure of the investigated meta- 
materials. The periodic structure becomes important when the 
wavelength of the electromagnetic wave is comparable with 
the lattice structure of the material^^. We proposed a more 
general description of the LHM, based on the concept of a 
periodic effective medium (PEM). This method enables us to 
distinguish between the resonant frequency dependence cor- 
responding to Eq. ^and effects of the periodicity of the struc- 
ture. We apply the PEM method for the analysis of numerical 
data obtained by the transfer matrix method (TMM). 

The paper is organized as follows. 

In Section ini we first explain basic ideas of the homoge- 
neous effective medium (HEM). Special attention is given to 
the correction of the phase of the EM wave at the interfaces, 
which is crucial for any retrieval procedure. 

In Section|ni]we define and analyze one dimensional peri- 
odic structures. The analyzed medium consists of thin slabs 
of homogeneous LH material separated by slabs of vacuum. 
We show that the approximation of such periodic medium 
by a homogeneous one give us effective parameters e and 
/i which possess unusual frequency dependences, similar to 
those observed when we approximate meta-materials by a ho- 
mogeneous medium. This proves that the periodicity of meta- 
material must be taken into consideration in the analysis of the 
effective parameters. 

The periodic effective medium is analyzed in two different 
formulations: continuous (Section [ill A> and lattice (Section 
IIIIB> . The latter is more relevant for the analysis of numeri- 
cal data since all known numerical algorithms use spacial dis- 
cretization. 



In Section llVl we analyze transmission data, observed from 
numerical simulations of periodic lattices of SRR, LHM and 
cut wires. We map these structures to periodic effective media 
which consist of homogeneous slabs separated by vacuum. In 
this formulation, e and fi of the homogeneous slabs are free 
from any modifications of the resonant behavior. To show the 
role of the periodicity of the meta-materials more clearly, we 
also analyzed a lattice of SRR in which we filled the gaps of 
the SRR by a dielectric with very strong dielectric permittiv- 
ity. This decreases the magnetic resonant frequency so that the 
wavelength of the incident EM wave is 25 times larger than 
the lattice period. We show that effective parameters again do 
not possess any deviations from resonant formula ([Q. 

A discussion of the applicability of various proposed mod- 
els to the analysis of transmission data is given in SectionlVl 
We discuss how the periodicity and anisotropy of the structure 
influence the transmission amplitudes and, subsequently, the 
effective parameters of the meta-materials. Final conclusions 
are given in SectionlVTl 



II. HOMOGENEOUS EFFECTIVE MEDIUM 

For the one-dimensional plain wave scattering problem at 
a homogeneous finite slab it is straightforward to obtain the 
scattering formulae. For the transfer matrices Tq for a single 
slice of vacuum and Tgiab for a single slice of homogeneous 
material with the thickness d we find in wave-representation"*"* 

with the elements 

a(d) — cos(qd) + - ( z + —] sm(qd) (3) 
2 V zj 

m = ^ - ^) ^"^(«'^)- ^4) 

In the continuum formulation and for normal incidence the 
momentum q inside the slab is related to the momentum k 
in the vacuum by the index of refraction n{k) = q/k, the 
impedance z is defined by z = ^{uj)k/q = q/{e{uj)k) for the 
TE and TM mode, respectively. Here, fJ,{uj) and e{uj) denote 
the frequency-dependent complex permeability and permittiv- 
ity of the homogeneous medium. On the lattice, ie. when we 
are going to compare with TMM simulation results, we have 
to take the modified dispersion relations 2 — 2 cos(fc) — oj^ = 
in the vacuum and 2 — 2cos((7) — ii{u!)e{uj) uj^ — in- 
side the slab into account. Then we have a modified q — 
acos (1 — /xe(l — cos k)) which gets noticeable at higher fre- 
quencies. Using the interrelation between the transfer matrix 
and the scattering matrix which defines the transmission (t±) 
and refection (r^) amplitudes. 
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we can calculate the transmission and reflection amplitudes 
for a sample composed of a left vacuum slice of length a, 
followed by homogeneous unit cells of length L in prop- 
agation direction, and terminated by a right vacuum slice of 
length b, 
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t_ = 
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In order to relate to the simulated scattering amplitudes com- 
puted numerically by the TMM by decomposition of the EM 
waves in the vacuum right of the sample with respect to the 
vacuum wave base left of the sample, it is convenient to in- 
troduce the normalized scattering amplitudes T and R which, 
after N unit cells, take the form 



T 



t 
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R = (3{-d) e-'^^"-^'^ T. 



(8) 
(9) 



In the continuum the scattering amplitudes of the homoge- 
neous slab are typically defined from interface to interface of 
the sample, ie. assuming a — h — Q. In the numerical sim- 
ulation this is not possible because of the lattice: we always 
have to make 1/2 vacuum-transfermatrix step from the last 
left vacuum slice into the sample and another 1/2 vacuum- 
transfermatrix step out of the sample onto the first right vac- 
uum slice. Therefore, the TMM scattering amplitudes T™"' 
and R '™"' are related to the normalized T and R involving an 
additional vacuum-phase compensation T = e^*'"' T'™"' and 
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Now we can resolve the above scattering formulae with 
given amplitudes T and R obtained from the simulation (or 
measurement) of a meta-material with respect to the mate- 
rial parameters impedance z{u!) and index of refraction 
If the solutions are (virtually) independent on the length of 
the sample those parameters define the homogeneous effec- 
tive medium (HEM) representation (or approximation) of the 
respective meta-material. Then we haveii 
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m, (11) 



with m G Z. Note that we obtain Zctt and rics from the scat- 
tering amplitudes only up to a common sign and the real part 
of the effective index of refraction, Rc Ucs only as a residue 
class. The former issue can be resolved by imposing addi- 
tional physical requirements, for instance Re z > (causal- 
ity). The problem of the residue class for Re rieff can be ad- 
dresses by considering different length {Li}, i G I C Z. 
Then we obtain a system of linear congruences, the solu- 
tion of which - if any - is a reduced residue class modulo 
2tt/ {k gcd {Li}) given by the greatest common divisor of the 
lengths Li. Since due to the inherent periodic structure of real 
meta-materials in simulations and experiments the lengths of 



the sample can only be integral multiples of the unit cell's 
length, the minimum possible ambiguity for Re ncfr will be 
a residue class modulo 271 /(kL) where L is the length of a 
single unit cell. For physical reasons we can assume a smooth 
frequency dependence between resonances which enables us 
to obtain Re noff(a;) as the corresponding residue class of 
piecewise continuous functions. The correct branch then has 
to be chosen exploiting additional physical information or as- 
sumptions of the model like the behavior of Uesi'^) at the 
plasma-frequency, in resonance induced transmission gaps an 
periodicity induced band gaps (discussed later). For known 
Ucfii'^) and Zcs{ljj) the effective permeability /i and permit- 
tivity e can be defined as 



eoff(w) = ncff(w)/2:off(w), 



(12) 
(13) 



respectively. 



Results for the effective parameters of the HEM approxi- 
mation of simulated meta-materials like arrays of SRR or cut- 
wires, LHM and even multi-gap SRR have been published by 
several authors ^i^iii'iSiSSi^iS. They all expose details which 
are in conflict with the simple effective medium behavior in 
term of a resonant and a plasmonic e(a;), originally pro- 
posed by Pendry, even under the assumption of an additional 
electric response of the SRR. Typical examples are also shown 
in Figs. |5] |6l |9] and ^| All results show resonance/anti- 
resonance coupling in jicsi'-^) and ecff(^) accompanied by 
negative imaginary parts, apparently different resonance fre- 
quencies for TT-cff and Zcff , the cut-off of the expected resonant 
positive (SRR) or negative (LHM) index of refraction, a mis- 
shapen, strongly asymmetric anticipated magnetic resonance 
in ^ for the SRR and LHM or electric resonance in e for the 
cut-wire, and finally a lot of unexplained additional structure 
(erratic stop-band and passbands) at higher frequency. 

Our extensive numerical simulations suggested that com- 
mon cause for all these problems has to be sought in the inher- 
ent periodicity, always present in the artificial meta-materials 
as they are composed of a repetitions of a single unit cell. To 
prove that the behavior is generic and really independent on 
the details of the unit cell, and that we can reproduce each of 
the effect above purely as a consequence of periodicity in the 
propagation direction, we investigated the most simple model 
for an effective medium with a non-trivial periodicity. 



III. PERIODIC EFFECTIVE MEDIUM 

To study the impact of the periodicity, or more precise the 
reduced translational symmetry of the sample in propagation 
direction, we consider a sample composed of a repetition of 
the unit cell shown in Fig. Q finite in direction of propaga- 
tion and infinite perpendicular to it. The unit cell consists of 
a thin homogeneous core of thickness d characterized by ar- 
bitrary ij.{lo) and £(a;), sandwiched by two slabs of vacuum 
with thickness a and h which break translational invariance. 
L is the length of one unit cell, N the number of unit cells 
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FIG. 1 : The layout of the single unit cell (a) and of a finite slab of 
the model periodic medium are shown. The shaded regions indicate 
the homogeneous core of the width d which is characterized by the 
chosen appropriately model functions and e(tj), sandwiched by 
two vacuum slabs. L is the length of a single unit cell, TV the number 
of unit cells in the slab in propagations direction. Periodic boundary 
conditions apply in the directions perpendicular to the propagation 
direction z. 



in propagation direction. To make a connection to our meta- 
materials we choose a simple Lorentz-type resonant form of 
fi{u!) and/or e(a;) to represent the magnetic and cut-wire re- 
sponse of the SRR. To model the LHM we would add a plas- 
monic term in e{uj) to account for the response of the continu- 
ous wires. Now we can calculate the scattering amplitudes for 
this model and subject them to the HEM inversion discussed 
in the previous section. The description (or approximation) of 
the scattering amplitudes for a given meta-material in terms 
of the effective parameters of such a periodic medium as de- 
fined in Fig. [2 will be denoted an "periodic effective medium" 
(PEM). 

The following results will show that this periodic medium 
can expose all the problematic effect discussed above. In a 
subsequent section we shall then demonstrate that this also 
applies to the simulated real meta-material. Their effective 
behavior can be decomposed into a "well-behaving" effec- 
tive response of the resonances and a contribution of periodic 
structure described by the PEM. 



A. Continuum formulation 

With the transfer matrices Tq and Tgiab introduced above, 
we can express the total transfer matrix of a finite slab of the 
periodic effective medium defined in Fig.^in the form 



= T^\NL) [To(6)T,i,b(d)To(a) 



N 



As expected from the z-inversion symmetry both transfer ma- 
trices To and Tgiab are unimodular, obviously is dot To = 1 
and a short calculation yields dctTgiab = a{d) a{—d) + 
0^{d) = 1. Therefore we can easily calculate the A^-th power 
of the unimodular 2x2 matrix above by diagonalizing it and 
computing the A^-th power of its eigenvalues^^. Using the in- 
teiTelation between the transfer matrix and the scattering ma- 
trix, we obtain the transmission and reflection amplitudes cor- 
responding to those computed numerically by the TMM, 



-ikNL 



(14) 
(15) 



Here, the Um = Un ip) are the Chebyshev polynomials of the 
second kind, Un{z) — sin((?i + 1) acos z)/(l — z^)^/^, taken 
at the argument 



P 



i{qd) cos [k{L — d)] 



(16) 



z + 



sm{qd) sin [k{L — c?)] 



The wave vector q = n{uj) k and the impedance z{u}) refer 
to the homogeneous core of the unit cell. For the normalized 
scattering amplitudes T and R after N unit cells we find 

T = (a(-d)e-^'=("+'') C/Ar_i(p) - UN-2{p)y' (17) 

R = l3{-d) e-''=("-*) Un-i{p) T. (18) 

Now we shall discuss what happens if we try to approxi- 
mate the explicitly periodic medium discussed above by a ho- 
mogeneous effective medium. This basically corresponds to 
our previous attempts to describe the periodic meta-materials 
by an homogeneous effective medium. We have two options: 
First, we could simply consider the analytic scattering am- 
plitudes (1171 I18> derived above to be those of a homoge- 
neous system and try to solve for effective material parame- 
ters EceiN, Lu) and /icff (A^, lu). This has the advantage that the 
approximation can deal with a possible residual length depen- 
dence of the approximate homogeneous medium, leaving an 
explicit possibility to assess the quality of the approximation. 
The disadvantage is that we have to handle the rather com- 
plicated structure of the formulae arising from the Chebyshev 
polynomials. The second approach is to assume that an ex- 
act correspondence of the periodic effective medium to an ho- 
mogeneous effective medium exists. This assumption is sup- 
ported by the length-independence (after appropriate phase 
compensation) of the conventionally inverted simulation data. 
If there is such a homogeneous effective medium we can write 
the transfer matrix of the periodic medium in terms of the 
transfer matrix for the homogeneous slab, 



T^\NL) [To(6) T,iab(d) To (a) 
- To-i(iVi)(^«)T,i,,(A^L), 



1 JV 



(19) 



which implies in particular for a system length of only a single 
unit cell 

T„(6) T,,,,(d) To(a) = (^^^^T^iabl^) • (20) 

Since for a homogeneous slab the identity T^^j^(i) = 
Tsiab(^-^) holds, finding a ^"^''T^^^^^L) that satisfies equa- 



tion i20\ in turn implies length independence of the homoge- 
neous effective medium description. Note that Tgiab (d) has 
only two independent elements, because f3{d) = —(3{—d) is 
anti-symmetric and the determinant a{d) a{—d) + 0^{d) = 
1 is fixed, such that we can calculate the matrix elements 
a(d) = (1 - 0^{-d))/a{-d) and (3{d) = -/3(-d) from 
a{—d), f3{—d). The assumption ( 12m imposes a restriction on 
the boundaries of the periodic medium in propagation direc- 
tion. The off-diagonal elements of Tgiab are anti-symmetric 
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FIG. 2: (color online) The HEM inversion fEQS. 1241 125> of the analytic continuum PEM scattering amplitudes fEQS. I17lll8> for model SRR- 
type material parameters uum ~ 0.13, LOmp ~ 0.16 for the magnetic, and LOe ~ 0.4, tOep ~ 0.5 for the electric response and 7 — 10"'* 
(cf. Eqs. 1291 130> . The homogeneous core located in the middle of the unit cell was d = L/10 thick. The retrieved real (red, purple) and 
imaginary (green, turquoise) parts of effective parameters are shown as a function of frequency uj. The dashed lines show the real (purple) and 
imaginary (turquoise) parts of the anticipated homogeneous parameters (Eos. 1291 130t and corresponding index of refraction and impedance. 
The dash-dotted black lines in Re rieff indicates the upper edge of the Brillouin zone, riedge = fcedge/A; = rmr/(kL). 



but on the left side of equation ( I20t this symmetry is broken 
by the phase factors g-ife(a-b) introduced in the 

off-diagonal elements by the two vacuum slabs. As a con- 
sequence the description as a homogeneous medium is only 
possible for a — h = Q. Besides choosing a symmetric unit 
cell in the first place we may alternatively compensate the 
factor Q-^Ha-b) -^^ j-j^g reflection amplitude R (which works 
simultaneously for all lengths, see equation [Tsl. effectively 
redefining the boundaries of the system such that the slab is 
centered in the unit cells. In terms of the (normalized) scatter- 
ing amplitudes T and R, for the single unit cell we then have 
the conditions 



f3{-d)T ^ R 



(21) 



(3,si-L) T. (22) 



We already know how to invert the right side of these equa- 
tions, this is just what we did in the retrieval procedure for 
the HEM in the previous section. Defining renormalized scat- 
tering amplitudes T' = t e-''''-^-'^^ and R' = Re-''''-^-'^\ 
we could apply the same procedure to the left side. Note that 
the possibility that we can always solve T and R for £cff(^^) 
and Hcsioj) guarantees a solution of equation j20t . In other 
words, there is always an exact, length independent descrip- 
tion of the periodic effective medium as a homogeneous ef- 



fective medium characterized by ecs{uj) and ^off ('j-')- There 
is no freedom to chose the boundaries of the homogeneous 
medium relative to the periodic medium. As shown above, 
we get the full information about the homogeneous effective 
medium which describes a given periodic effective medium 
characterized by n{k), z{k) and the geometry d, L already 
from the first unit cell. Inserting the renormalized transmis- 
sion and reflection amplitudes (1211 122> for a single unit cell 
into the inverted scattering formulae above we obtain 



cos (ncfr(fc) kL) 



i{nkd) cos [k(L - d)] (23) 
in(nfcc?) sin [k{L — d)] 



where n(k) and z{k) are the parameters of the material slab 
in the middle of the unit cell of the periodic effective medium. 
With the p defined in equation il6\ and q = n{k)k we can 
write simpler 



ncs(k) = ±— arccos 
kL 



p{n, z; k) 



27r 



with TO e Z. The problem with the signs of ricft and z^s, 
as well as with the ambiguity of Rc ?icff is similar, and can be 
resolved the same way as for the case of the homogeneous slab 
discussed above. Analogously we can express the impedance 
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Zcff of the effective homogeneous medium in terms of the n 
and z of the homogeneous core as 



Zeff{k) = ±\ 



l2p+ + {z- l/z)sm{qd) 



(25) 



(26) 



2p+ - (z - l/z)sm{qd) 
where q — n{k)k and 

= cos{qd) sin [k{L — d)] 

+ ^ + i j sin((;d) cos [k{L — d)] . 



The parameters of an effective homogeneous medium de- 
scribing the periodic material from Fig. [2 which have been 
obtained from the formulae (I24> and (I25> . are shown in Fig.|2] 
for a concrete example of SRR-type pi{oj) and £(w). For the 
homogeneous core in the middle of the unit cell we have cho- 
sen 



e{k) = 
Kk) = 

with model functions 



l + (i/d)[£ff(fc)-l]. 



and 



eH(cj) = 1 - 



■ 



(27) 
(28) 



(29) 



(30) 



to emulate the anticipated magnetic and electric resonances 
of the SRR. For a LHM-type behavior we have to add the 
plasmonic response of the continuous wire in the permittivity, 



origination from the periodicity rather than from the under- 
lying material properties. The qualitative behavior presented 
in Fig.|2is generic for a wide range of parameters (jj,„, LOmp, 
ujf,, ... of the resonances and L, d of the geometry. If the peri- 
odic medium model is used with only the electric resonance or 
with an additional plasmonic term in e(uj\ it qualitatively re- 
produces the observed deviations from the expected plain ef- 
fective medium behavior published for the array of cut-wires 
and the LHM, respectively. 

Although the curves show most of the discussed abnormal- 
ities in the HEM description of the SRR, the analytic descrip- 
tion matches the simulation and inversion results for the real 
meta-material present in literature not in all aspects. Clearly, 
there are problems very close to the resonance frequencies. 
Instead of the divergence in the effective index of refraction 
rioff being virtually cut off at the upper edge of the first Bril- 
louin zone as observed in the simulations of the actual SRR 
meta-material, the analytic description produces a series of 
consecutive band gaps at the boundaries of the first and higher 
Brillouin zones and a lot of structure in the imaginary part 
of ricff . The same holds for the analytic description applied to 
the periodic effective medium model of the LHM (not shown). 
Here, we particularly miss the cutoff at the bottom of the neg- 
ative ricft region. In either case the underlying lattice in the 
simulation starts to become visible. Since the lattice has a 
finite lattice constant it cannot support arbitrarily large mo- 
menta, such that we expect additional effects if the continuum 
momentum q reaches the order of vr/aiatticc- In order to un- 
derstand also the details of the retrieved HEM parameters in 
our simulation of real SRR and LHM meta-materials we have 
to take the discretization lattice of the employed TMM into 
consideration. To see the modification of the continuum re- 
sults by the discretization lattice we have to derive the scatter- 
ing formulae for the periodic medium model on the lattice. 



(31) 



According to a simple effective medium picture, we would ex- 
pect that we can approximate a homogeneous unit cell char- 
acterized by /iff (til) and eni'-^) by concentration the mag- 
netic and electric polarizations into the homogeneous core 
of the periodic medium. Fig. |5] shows the actual effective 
impedance z^si'^) and index of refraction nog(a;) obtained 
via the HEM inversion of the periodic medium. Comparing 
with the expected effective medium behavior (dashed lines) 
we clearly see the typical anomalies in the shape and positions 
of the resonances, the same qualitative behavior as observed 
for real SRR meta-materials in the literature and our own pre- 
vious work. The effective parameters of the HEM approxima- 
tion of our periodic medium model show the resonance/anti- 
resonance coupling in /ioff (tij) and eoffCf-^) together with the 
negative imaginary part Im /ioff(w) < around the mag- 
netic resonance frequency cj,,,, and also a very involved be- 
havior close to the cut-wire resonance lo^. The effective in- 
dex of refraction is cut-off at the edge of the Brillouin zone 
which corresponds to the appearance of additional band gaps 



B. Lattice formulation 



We follow the TMM introduced for the Maxwell equations 
by Pendry^'i^ii^SiS in the formulation described by Markos 
and Soukouli&^. The electric and magnetic field, together 
with the spatially dependent material relative constants /^rci (r) 
and £rci(r) which define the meta-material, are discretized 
on the bonds of mutually dual lattices {m} and {rii}. With 
the renormalized material constants ei(m) ~ iuj Eq erci(ni + 
ei/2) and /ii(m) = ioj fiQ /ii.ci(m — ej/2), used throughout 
this section, we can write the transfer matrix equations for a 
stratification in z-direction for the two independent compo- 
nents i G {x, y} of the electromagnetic field. Using (quasi- 
)periodic boundary conditions in the _Lz plane we can intro- 
duce a Fourier representation of the fields with respect to this 
plane defining an in-plane momentum q. To derive a scatter- 
ing formula corresponding to the continuum case considered 
in the previous section we restrict ourself to the most simple 
case of normal incidence, ie. zero in-plane momentum q = 0. 
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Then the transfer matrix for normal incidence takes the form 
H ) = (32) 

B{m, + 1) l + B(m^ + l)A(m,) J [^U 

The generally q-dependent matrices A and B reduce to sim- 
ple off-diagonal form, with the product B(m2 + l)A(m2) 
diagonal, 



B(m,) 



IJ-yirriz) 
-IJ-ximz) 

-Sylmz) 
Sxirriz) 



(33) 
(34) 



such that the transfermatrix i32\ factorizes, reorder- 
ing the electromagnetic field vector in the form 
{Ex, Hy, Ey, ~Hx)^ , into a two-fold degenerated block- 
diagonal structure 



Ex Ey 

1 fJ^yirriz) 
Sxiniz + l) 1 + ex{mz + l)fiy(mz) J \ Hy 



(35) 



Ex Ey 



Without loss of generality we can restrict ourself to consider 
just the first polarization. We denote the single-polarization 
transfer matrix for the q = modes in the last equation 
T{mz). It is expedient to introduce the decomposition 

T{mz) = T^iniz + 1) Tf_,{mz) (36) 
1 ^ 1 fi{mz) , ^^^^ 



e{mz + 1) 1 / I 1 



Further we can factorize the r into a vacuum and a mate- 
rial contribution, = '''(e-£vac)''£vttc — ''evac'''(e-evac) 
related to the polarization for the magnetic and analog for 
the electric field step. Note the renormalized vacuum per- 
mittivity Evac = iujEo and permeability /ivac = ioJl^o- As 
expected, T{mz) is unimodular Now we can easily find the 
eigensystem; the eigensystem of the vacuum transfer matrix 
defines the plain wave basis on the lattice which we use to 
define the scattering formalism. Because of the unimodu- 



larity the two eigenvalues A 



are mutually reciprocal 



and for the propagating modes we are interested in on the 
unit circle, ie. k is real. We get the characteristic polyno- 
mial — A [2 + ex{mz + 1) ^y{mz)\ + 1, hence cos k = 
1 + £x{mz + 1) iJ,y{mz)/2. The two signs of k correspond to 
the right- and left-moving waves. Note that e and /i implicitly 
contain the oj dependence. To obtain the scattering matrix on 
the lattice we need the wave-representation of the total trans- 
fer matrix of a unit cell. The right and left eigenvectors of 

T{mz) are distinct, R\{mz) — (l, (A — 1)/ ^jLy['mz)) and 
Lx{mz) = (1, (A - l)/ex{mz + 1))*^/(A + 1)*, and satisfy 
the orthogonality relation L^, {rriz) Rx- (mz) ~ Note that 
we applied the common normalization to the left-eigenvectors 



in order to normalize the electric field component of all right- 
eigenvectors to one. This is requked for a clean definition of 
the scattering amplitudes analog to the continuum case. Fur- 
ther, the two right- and the two left-eigenvectors are linearly 
independent, respectively. Therefore we may group the two 
right- and the two left-eigenvectors of the vacuum transfer ma- 
trix into the matrices 



/ 1 

V ' 

Rn = 



1 e-''' - 1 



1 



tUJ£0 
— ik ^ 



(38) 



(39) 



where the eigenvalues A ~ satisfy the vacuum dispersion 
relation 2 — 2 cos k + /ivac£vac = 2 — 2 cos k — lu'^hoEo = 
for the vacuum wave vector k, and use the projector Rq 
to obtain the wave-representation of the total transfer matrix 
Ttot of the finite system as 



Ttot(fc) = L+(fc)TtotRo(fc)- 



(40) 



Then we get the usual definition of the scattering amplitudes 
from the correspondence between the scattering and the trans- 
fer matrix given by equation Q. 

The homogeneous slab. Now we have to consider the total 
transfer matrix of our meta-materials. The most simple case 
is just a homogeneous slab of finite length. On the lattice, the 
composition of the total transfer matrix depends on the ma- 
terial discretization. We compute the total transfer matrix by 
starting from a right-eigenvector of the vacuum base at the last 
vacuum site = just before one side of the sample and ap- 
ply successively the single-step transfer matrices T{mz) until 
we reach the first site niz = n + 1 right of the sample for 
which the T{mz) is a vacuum step again. We have n ma- 
terial layers inside the sample but n + 1 transfer matrix step 
which depend on the material parameters ^y, Ex of the sample. 
Since we only have to consider a single polarization, we drop 
in the following the y, x-indices in ^ and e in order to improve 
readability. Because in the discretized Maxwell equations the 
electric and magnetic fields live on mutual dual lattices, we 
distinguish three different single step T{mz) inside the sam- 
ple instead of only one, as one would expect for a homoge- 
neous slab. T{mz) depends on fJ,{mz) and E{mz + 1). There- 
fore the first step r(0) inside the sample sees only the electric 
response £(1) but no magnetic response of the material. The 
subsequent steps see both, e and /i, and are constant across the 
bulk of the sample. The last step back into the vacuum behind 
the slab is special again. Both steps across the boundaries 
of the sample depend on the chosen material discretization. 
Here we adopt a symmetric material discretization^^ which 
respects the z-isotropy such that the the steps into and out off 
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FIG. 3: (color online) The HEM inversion (Has. 1101 II it of the analytic lattice PEM scattering amplitudes (Has. I46> for model SRR-type 
material parameters ujrn = 0.13, LUmp = 0.16 for the magnetic, and uje = 0.4, cuep = 0.5 for the electric response and7 = 10"'' (cf. Eqs.|^ 
l30l . The homo geneous core located in the middle of the unit cell was d = L/10 thick. The retrieved real (red, purple) and imaginary (green, 
turquoise) parts of effective parameters are shown as a function of frequency lo. Note the reduction of the multiple band gaps seen in Fig.|2| 
around the the resonances to a single gap before each resonance. The dashed lines show the real (purple) and imaginary (turquoise) parts of 
the anticipated homogeneous parameters (Has. 1291 130> and corresponding index of refraction and impedance. The dash-dotted black lines in 
Re Weff indicates the upper edge of the Brillouin zone, n^dgc = fccdgc/fc ~ mn /(kV). 



the sample become equal. Then we may calculate the wave 
representation as 



= L+ T^X_ R [ L+T,T^ R ] " L+ T,^,_ Tvac Ro 



where {L^ ^ /Jq} is the eigenbase of the vacuum transfer ma- 
trix step Tvacuum with the eigenvalues Aq = e^**^ as before, 
but {L"'", R] now denotes the eigenbase of the transfer matrix 
step TgT^ inside the homogeneous medium with the eigenval- 
ues A = e±'«. We made use of the afore mentioned identity 
Te^Te^ = Tj^+e^. The Symmetric material discretization intro- 
duces the averaged e = (e + £vac)/2 at the material's surface. 
As shown above, the wave vector in the vacuum k and in- 
side the homogeneous slab q satisfy the dispersion relations 
2 — 2 cos k + /ivac^vac — and 2 — 2 cos q + fie = 0. Since 
the matrix L+r^r^R in Eq.|^is diagonal, we basically have 



to calculate the matrix L+r^ 
obtain for the homogeneous slab 



^Rq. After some algebra we 



Ttot(fc) 



^corc ''^) 



vacuum 



(k) 



with the diagonal Tvacuum (fc) and 



(42) 



1 xnG{X)-^G{{) (A"-3i,)(2-G(l)) 
C V -(^" - A^) (2 - G(l)) -A" G(i) + ^ G(A) 



where 



G(A) 
a{k) 

c 



A 

Ao(fc) 



Xa{k) -i 

(A-A-i)(Ao-Ao^) ^ 

^Mvac 



— fie a{k) a{~k), (43) 
— , (44) 
(45) 



with a{k) + a(— fc) = 1 and consequently G(l) = 1 — 
jie a{k) a{~k). Further we have a{k) — a{—k) — C/(A — 
A"^) and G(A) = G(A~^) + C- Note the anti-symmetry of 
the off-diagonal elements. Using again the definition of the 
scattering matrix (|5}, we find the transmission and reflection 
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FIG. 4: (color online) The HEM inversion ("Eqs. IIOIIIH of the analytic lattice PEM scattering amplitudes (Egs. I46> for model LHM-type 
material parameters Um, = 0.13, lOm-p = 0.16 for the magnetic, and lo^ — 0.4, lOep ~ 0.5, ojp — 0.27 for the electric response, and 7 = 10^* 
(cf. Eqs. 1291 13 1> . The homogeneous core located in the middle of the unit cell was d = L/10 thick. The retrieved real (red, purple) and 
imaginary (green, turquoise) parts of effective parameters are shown as a function of frequency uj. The dashed lines show the real (purple) and 
imaginary (turquoise) parts of the anticipated homogeneous parameters (Eqs. 1291 131> and corresponding index of refraction and impedance. 
The dash-dotted black lines in Re riefi indicate the edges of the Brillouin zone, ncdgc = fcodgc/'^ = mn /(kL). 



amplitudes as 



Te^+£^ T^^Te^ and consequently Tg-e Tg-e^^, = 1 



CAo 



(46) 



A-"G(A)- A"G(A-i) 
r+(c^,n) - (CAo)-MA"-A-")(2-G(l))t_(a;,n). 

The non-vacuum factor of the lattice transfer matrix (I42t ap- 
pears to have the same symmetries as the transfermatrix of the 
homogeneous slab in the continuum: the off-diagonal terms 
are anti-symmetric, the diagonal terms are mutual complex 
conjugates if £rci and /iici are real. 

The periodic medium. Knowing the transfer matrix of the 
finite slab it is now easy to obtain the transfer matrix for a sam- 
ple of multiple unit cells of the homogeneous as well as the 
periodic medium with the unit cell corresponding to Fig.^. 
We can reduce the wave representation of the total transfer 
matrix to a product involving the wave representation of the 
homogeneous core we already know and some additional vac- 
uum transfer matrix steps Tq for the free space in the unit cell. 
We assume the measures a, h and d in Fig.^to correspond to 
fia, rih and n^; layers on the lattice. Then we get for the total 
transfer matrix of N unit cells of the periodic medium using 



(47) 



N 



N 



with the Tcorc {k, Ud) defined in equation \A2\ . Since the phase 



factors A, 



-(n5-n„) 



and Aj;"-"" 



introduced by the two vacuum 



slabs in the bracket on the last line of equation ( I47> do explic- 
itly break the anti-symmetry of the off-diagonal elements that 
is present for the single homogeneous slab in the continuum 
and, in the symmetric material discretization, also on the lat- 
tice, we can obtain a representation of the periodic medium 
by a homogeneous medium only for the case = Ua- As 
already explained for the continuum case this is not a real re- 
striction but instead just fixes the definition of the effective 
boundaries of the periodic medium. In the numeric simulation 
we have to explicitly compensate the corresponding vacuum 
phases in the scattering amplitudes. We can use the Cheby- 
shev formula to explicitly calculate the iV-th power such that 
we get the transmission and reflection amplitudes for the pe- 
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riodic medium after N unit cells in propagation direction as 



An 



CK 



nb+ria 



Un- 



Un- 



r+{u;,N) = 

^Xn. _ A-"-) (2 - G(l)) f/jv-i t^{uj,N) 



(48) 



(49) 



CA- 



An 



where the argument of the Chebyshev polynomials Un (p) is 
given by 



P 



cos{qnd) cos [k{na + Ub)] 
G(A)+G(A-i) 



(50) 



c 



siii{qnd) sin [k{na + Ub)] . 



As for the continuum formulation, we actually get all the in- 
formation about the meta-material from the single unit cell. 
Comparing the scattering amplitudes j48l I49t on the lattice 
with the normalized scattering amplitudes for homogeneous 
slab in the continuum tells us how to do the phase compen- 
sation for the lattice-tmm results: T = Xq^ t- and R = 
The condition for T arises from the additional 



A, 





vacuum step Tq into the slab on the lattice, the compensation 
in R results from the symmetric definition of the boundary of 
the unit cell which is required to describe the periodic by a 
homogeneous medium as explained above. 

Continuum HEM inversion. Again we ask whether the 
model periodic medium from Fig.^can be represented by an 
effective homogeneous medium. Here we have two choices: 
{i) we can compare the scattering amplitudes of the lattice pe- 
riodic medium with the scattering formulae ( 1211 \22\ derived 
for the homogeneous slab in the continuum, or we can in) 
compare with the lattice scattering formulae for the homo- 
geneous slab derived in this section. Moreover, we have to 
decide which material discretization to use. In this paper we 
will concentrate on comparing the lattice scattering results to 
the continuum scattering formulae for the homogeneous slab, 
as we previously did with the standard inversion procedure 
to obtain effective ecsi'^) and /ioff(w) from the meta-material 
simulations. 

Analytically, the effective material parameters obtained 
from the HEM inversion for the lattice formulation of the 
model periodic media used in the last section to emulate the 
SRR and LHM meta-materials are shown in Fig.|3]and Fig.0] 
respectively. As expected, the qualitative behavior is very sim- 
ilar to that found with the continuum formulation. All the 
problematic effects seen in the previously published simula- 
tions, like resonance/anti-resonance coupling, negative imag- 
inary parts, deformed resonances, bad gaps and so on, are 
present. The major difference to the continuum formulation 
becomes visible around the resonances. Where we previously 
found a series of tiny periodicity band gap around the res- 
onances, in the lattice formulation we obtain a much sim- 
pler structure with basically one gap before each resonance. 
This is in excellent agreement with the numerical simulations. 



HEM Homogeneous Effective Medium, a homogeneous 

medium characterized by ^i{i-o) and e{ijj) which, 
substituted for a finite meta-material slab, length- 
independently reproduces (or approximates) the 
given scattering amplitudes. Here always used in 
continuum formulation. Finding a HEM for given 
T, R is called HEM inversion (if exact) or HEM 
approximation. 

PEM Periodic Effective Medium, a most simple peri- 

odic model-medium defined by /icoro(<^), ecoro(<j-') 
and a geometry shown in Fig. Q which length- 
independently reproduces (or approximates) given 
scattering amplitudes. Also used a priori with given 

/-^corc t core 

(oj) to demonstrate effects of the pe- 
riodicity. Here used in lattice formulation. 
HEM(PEM) The HEM which reproduces the scattering ampli- 
tudes calculated analytically from a given PEM. 

TABLE I: Summary of the effective medium related acronyms. 



hence, expectedly, the lattice formulation compares much bet- 
ter to numerical simulations also obtained via discretization of 
the Maxwell equations than the continuum formulation. The 
discussion of further details we shall defer to a dedicated sec- 
tion below. 



IV. SIMULATION RESULTS 

In this section we now present actual TMM simulation re- 
sults for real SRR and off-plane LHM meta-materials. All 
numerical simulation are done using an implementation of 
the TMM method described by Markos and Soukoulis'^. The 
meta-materials are uniformly discretized on a cubic lattice us- 
ing a symmetric material discretization. The dimensions of 
the unit cell are 6x10x10 mesh steps, the single-ring SRR is 
a square ring of 7x7 mesh steps with a gap in the top side 
one mesh step wide. Propagation is for all cases along the 
SRR plane with the polarization of the incident plain wave 
such that the electric field is parallel to the two continuous 
sides of the SRR. Therefore we have only magnetic coupling 
to the magnetic resonance of the SRR^"''*. Periodic boundary 
conditions apply to both directions perpendicular to the di- 
rection of propagation. For the off-plane LHM we add a one 
mesh-step thick continuous wire in front of the SRR such that 
the position of the wire is symmetric in the middle between 
two periodic repetitions of the SRR plane and centered with 
respect to the gap in the SRR. The direction of the wire is par- 
allel to the continuous sides of the SRR, thus parallel to the 
incident electric field. All components of the meta-materials, 
the ring ring of the SRR and the continuous wire, are made 
from metal characterized by a constant relative permittivity of 
Emctai = (-3.0 + 5.88 i) 10^ and /imctai = 1- Note that the 
results do not depend much on £mctai as long as it does not 
fall below a certain threshold'^. The chosen value is reason- 
able to emulate metals like Cu, Ag, Au in the range of GHz 
to a few THz. The rest of the unit cell is vacuum, there are 
no dielectric boards. The special geometry of the unit cell has 
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been carefully chosen to preserve the inversion symmetry of 
the unit cell in the two directions perpendicular to the direc- 
tion of propagation. This allows us to consider the scattering 
for only one polarization as it avoids complications by cross- 
polarization terms in the scattering amplitudes'^. In this pa- 
per, we concentrate our consideration on the region around the 
magnetic resonance frequency curn, where we expect /icg (w) 
to become transitionally negative, for two reasons: first this is 
the region of interest for any left-handed application, and sec- 
ond, this is the frequency window for which simulation data 
is typically shown in the literature. A more detailed investi- 
gation of the higher frequency region, particularly the vicinity 
of the of the electric cut-wire response of the SRR and the in- 
termediate periodicity band gaps will be published elsewhere. 

In the following we show HEM inversion results for the 
scattering data numerically obtained for the meta-material 
with the TMM. After the correct vacuum-phase compensation 
described above the inverted HEM scattering formulae dlOl 
II 1> are applied to the simulated T and R for meta-material 
slabs with a thickness of one, two and three unit cell in prop- 
agation direction. We shall denote the results as nhem(w) and 
2;hcm(^^) or fihcm = '^hem 2^hem and Ehom = '^hem/^^hem, Cor- 
respondingly. This approach is the same as chosen in the lit- 
erature. Then we find the PEM approximation for the simu- 
lated meta-material using the lattice formulation for the ana- 
lytic scattering formula of a model periodic medium consist- 
ing of a homogeneous core which is a single discretization 
mesh-step thick and located the unit cell in the plain of the 
SRR gaps and the LHM continuous wire. This constitutes 
the lattice equivalent of a single scattering plain in the con- 
tinuum. A model periodic medium, characterized by effective 
material constants and eito) of the homogeneous core, 
which reproduces the simulated T and R independent on the 
system length is called a periodic effective medium. The nu- 
meric inversion of the lattice scattering formulae ( 1481 149> is 
applied to the simulated T and R for the firsts unit cell of 
each meta-material, providing us with effective material con- 
stants /icorc(w) and ecoro('^) for the homogeneous core of the 
PEM approximation. From the core parameters we can derive 
two further sets of effective parameters. First, we calculate 
the HEM inversion of the PEM scattering data obtained from 
the retrieved /icore('^) and £core(w) and compare the results 
with the HEM inversion of the direct simulation data to as- 
sess the quality of the PEM approximation. We denote this 
as HEM(PEM). Second, we introduce the material parameters 
/^pcm('^) and epom(^) of a homogeneous unit cell that would 
correspond to the PEM approximation in the effective medium 
limit, equating the total electric and magnetic polarizations of 
the respective unit cells, 



Mpcin('^) 



1 



na + nd + Ub 



and 



epcm(t^) = 1 + 



nd 



na + nd + lib 



(^corc(t^) ~ 1) (51) 



(ecorcM-1). (52) 



The idea of this definition is to obtain parameters which we 
can compare with those of the HEM inversion, becoming 



equivalent with the latter if we can truly neglect the period- 
icity of the material. This allows us, to some degree, to con- 
sider the meta-material' s electromagnetic response as being 
composed of an actual contribution of the internal geometry 
of the meta-materials constituents and an exphcit contribution 
of the periodic arrangement. 



A. SRR 

From the naive effective medium picture we expect the SRR 
to expose at the magnetic resonance brought about by the LC- 
oscillator type response of the split ring to an perpendicu- 
lar external magnetic field a Lorentz-type resonant form/ in 
the permeability iics{^) but a unaffected, weakly frequency 
dependent permittivity ecs{oj). If the magnetic resonance is 
strong enough, ie. the imaginary part in fi is small, we should 
find a isolated region where Re /Xoff (w) < 0. The HEM ap- 
proximation of the actual simulation data is shown as effective 
impedance z\icm{^) and index of refraction nhcm('^) in Fig.|5l 
as the colored points. The different colors (orange, yellow, 
green) correspond to different lengths of the meta-material of 
one, two and three unit cells in propagation direction. We see 
the typical behavior around u!,n we are already familiar with 
from previous work^'. Instead of the expected form, we ob- 
tain a sharp cut-off in Re nhem('jj) on the low frequency side, 
accompanied by a non-zero Im nhcm(w). The adjacent region 
with Re rihom ~ and significant Im nhom > should corre- 
spond to the negative produced by the magnetic resonance. 
Even more disturbing is the fact that the peaks in Zhcm('jj) 
which should coincide with the peaks in rihcm('-^) directly at 
the resonance frequency ojm do appear at substantially lower 
frequency. This behavior makes it hard to determine ijj,„ for 
the meta-material since the retrieved nhcmlt^) and Zhem('jj) 
mutually disagree upon the value. For the first time we show 
that the HEM approximation of the meta-material is, apart 
from some additional noise for longer systems, indeed length 
independent. Although we only show data for the first three 
unit cells we confirmed the length independence for all system 
lengths up to 11 unit cells. The dash-dotted line in Fig. [Sja) 
indicates the upper edge of the first Brillouin zone, ricdge — 
fcodgo/fc — Tr/{kL) where L is the length of the unit cell in 
propagation direction. Clearly the cut-off of Re n\icm{^) co- 
incides with this line. Note also that the peaks in Zhcm('jj) 
appear exactly when Re riiicm('^) reaches ncdgc(^)- This be- 
havior is generic, qualitatively the same is observed for dif- 
ferent sizes of the unit cell and different geometries of the 
SRR, including single-ring and double-ring SRRs as well as 
more symmetric multi-gap SRRs"^. The corresponding effec- 
tive permittivity and permeability of the HEM approximation 
are shown in Fig. |6l Again, the colored points represent the 
simulation data for one to three unit cells. The most strik- 
ing deviation from the expected effective medium behavior is 
the resonance/anti-resonance coupling between Re //hem(w) 
and Re Shcmi'^), accompanied by a significant negative imag- 
inary part of the permittivity, Im ehcm(^) < 0. Moreover, 
the negative region of the /i resonance is strongly but charac- 
teristically deformed and not ascending monotonically from 
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FIG. 5: (color online) For the simulated SRR meta-material the effective index of refraction noff(a;) and impedance Zc:ft{uj) are shown. 
The colored curves represent the HEM approximation (Egs. llOlllH of the simulation data for the first 3 unit cells, the solid black line the 
HEM(PEM) approximation, and the dashed line the PEM approximation (cf. Eqs. 1511 152> obtained from the first unit cell data. Note the 
different positions of the resonance for ncB(io) and Zob{uj). 





FIG. 6: (color online) For the simulated SRR meta-material the effective permittivity ecff(<j-') and permeability ij,ch{ui) are shown. The colored 
curves represent the HEM approximation (cf. Egs. llOllTTl of the simulation data for the first 3 unit cells, the solid black line the HEM(PEM) 
approximation, and the dashed line the PEM approximation (Eqs. I51||521 obtained from the first unit cell data. Note the anti-resonant behavior 
of the permittivity and the misshapen magnetic resonance in the frequency interval where nofi(tj) reaches the edge of the Brillouin zone. 
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FIG. 7: (color online) For the simulated low-frequency SRR meta-material the effective index of refraction Wcff (oj) and impedance (oj) are 
shown. The colored curves represent the HEM approximation (Egs. llOlll H of the simulation data for the first 3 unit cells, the solid black line 
the HEM(PEM) approximation, and the dashed line the PEM approximation (cf. Egs. lSllB^ obtained from the first unit cell data. Note that 
far away from the edge of the Brillouin zone HEM and PEM approximation, and the expected homogeneous medium behavior coincide. 
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FIG. 8: (color online) For the simulated low-frequency SRR meta-material the effective permittivity £eff(<^) and permeability /^off(<jj) are 
shown. The colored curves represent the HEM approximation (cf. Egs. llOlim of the simulation data for the first 3 unit cells, the solid black 
line the HEM(PEM) approximation, and the dashed line the PEM approximation (Eqs. l51ll52t obtained from the first unit cell data. In the 
low-frequency limit the resonance/anti-resonance coupling as well as the negative imaginary parts disappear. 
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FIG. 9: (color online) For the simulated off-plane LHM meta-material the effective index of refraction rieff (a;) and impedance Zcff (i^) are 
shown. The colored curves represent the HEM approximation (Egs. llOlll H of the simulation data for the first 3 unit cells, the solid black line 
the HEM(PEM) approximation, and the dashed line the PEM approximation (cf. Eqs. I51||521 obtained from the first unit cell data. 





FIG. 10: (color online) For the simulated off-plane LHM meta-material the effective permittivity ScS (i^) and permeability h^b (uj) are shown. 
The colored curves represent the HEM approximation (cf. Eqs. llOlfm of the simulation data for the first 3 unit cells, the solid black line the 
HEM(PEM) approximation, and the dashed line the PEM approximation (EQS. I51l|52t obtained from the first unit cell data. 
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a negative divergency. Of course, the divergencies of the ef- 
fective parameters would be blurred in the presence of large 
imaginary parts in /i or e, but at least for our simulations using 
almost perfect metals in vacuum we would expect reasonably 
sharp divergencies. As a consequence of the length indepen- 
dence of nhcm and Zhcm also the retrieved Re /ihem('jj) and 
Re ehcm(w) are basically length independent. 

In the analytic sections above we demonstrated that the pe- 
riodicity can produce all that kind of effects violating the ef- 
fective medium picture in our simulations. Now we show 
that the PEM approximation of the simulation data yields rea- 
sonable effective parameters free of the above artifacts. The 
dashed curves in Figs.|5]and|6lrepresent the effective param- 
eters fipcai{i^) and epom(^), with the corresponding Upcmi^) 
and Zpoin(i^), extracted from the PEM approximation of the 
simulation data as described above. For the PEM unit cell 
we chose ria — 5, nj, = 1 and ni, = 4 to fit the symmetry 
of the meta-materials unit cell that embeds the 1x7x7 SRR 
into a 6x10x10 unit cell. From the corresponding /icorc('-^) 
and ecorc(w) we can in turn calculate the scattering ampli- 
tudes of the PEM and subject those to the HEM inversion. 
The resulting HEM(PEM) parameters are shown as the solid 
black lines which virtually coincide with the first unit cells 
simulated HEM parameters, hence basically also with those 
of the the longer systems, proving that the PEM approxima- 
tion is reliable. Obviously the PEM parameters behave exactly 
as we expected from the effective medium picture in the first 
place. There is a clean resonance in npom(^) and the peaks of 
JT-pcml^) and Zpcm(w) now appear at the same frequency, ren- 
dering ujrn well-defined. There is no resonance/anti-resonance 
coupling between fJ.pcm{ijj) and epem(w) anymore. The per- 
mittivity is roughly constant across u!„i and does not show any 
negative imaginary part. The permeability /ipom(^) exposes a 
clear, anti-symmetric resonance at a;,„ in its real part in con- 
junction with a symmetric absorption peak in the imaginary 
part at the same frequency. Note further that the frequency 
i^mp where Re /Lthem(^) get back to positive values after the 
magnetic resonance is substantially smaller than the corre- 
sponding frequency for Re /ipom(w)- The procedure yields 
analogue results for different SRR geometries and different 
size of the unit cell (not shown). 

Effective medium behavior was expected in the first place 
because the vacuum wavelength of the incident radiation at 
the magnetic resonance frequency a;,„ is large compared to 
the size of the unit cell for the customary LH and related 
meta-materials. In most experiments and simulations this ra- 
tio is in the range of ten to five5''' 7-y-iO-i3.22.23,2736.37.38 ^^le 

SRR meta-material presented above we went intensionally to 
the lower extreme to make the effects of the periodicity more 
pronounced and more easy to identify. In Figs. and Owe 
present the analog HEM and PEM approximation results for 
a SRR with five times lower u!,n- The size of the SRR and 
the unit cell are kept constant. To lower the resonance fre- 
quency of a SRR without changing the size one would usually 
decrease the width of the gap in the SRR increasing its capac- 
itance. Due to the limitations of the used TMM implementa- 
tion (uniform discretization) this was not feasible. Therefore 
we adopted the alternative possibility to place some high di- 



electric constant material inside the SRR gap, which serves 
the same purpose and can be used to emulate a narrower gap. 
For this low-frequency SRR the wavelength to unit cell size 
ratio around uj„i is roughly 25, higher than in any published 
LHM. The behavior of the retrieved effective nhcm('jj) and 
^hcm('jj) in Fig. is now qualitatively as expected from the 
effective medium picture, the refractive index and impedance 
of the HEM approximation virtually coincide with the PEM 
parameters npein(w) and Zpeni(w). As the magnetic reso- 
nance is now shifted far below the edge of the Brillouin zone 
there is no visible cut-off in Re nhcm('jj) and the resonance 
peaks in n and z appear at the same frequency cj,,, . Note that 
/^hcm('jj) reaches unity away from u!,n to either side. Also 
the the effective permeabilities of HEM and PEM in Fig. |S] 
do coincide, exposing a clean anti-symmetric resonance in 
Re ^ and a symmetric positive absorption peak in Im /i. Sur- 
prisingly, though weak there is still a noticeable residue of 
the resonance/anti-resonance coupling left in Re £hcm('-^) to- 
gether with the corresponding negative imaginary part. Only 
here the HEM approximation deviates from the PEM approx- 
imation which expectedly does not show these effects. All 
effective parameters are almost perfectly length independent. 
This has been verified for up to 10 unit cells in propagation 
direction (not shown). Note also the absence of the additional 
noise observed in the effective parameters for longer systems. 
These results show that the artifacts in the HEM approxima- 
tion which we identified as effects of the meta-material' s pe- 
riodicity vanish if we approach the effective medium limit. 
At low frequencies HEM and PEM approximation converge, 
however, even in this extreme low-frequency limit remains of 
the periodicity artifacts are still visible. Since for virtually all 
the meta-materials measured or simulated that have been pub- 
lished the wavelength to unit cell size ratio is worse than 25, 
we expect the meta-material' s periodicity required to be taken 
into account. 



B. LHM 

The most simple way to obtain a left-handed material is to 
add an appropriately dimensioned^^'^'' continuous wire to the 
SRR considered above. From the naive effective medium pic- 
ture we expect for the LHM a Lorentz-type resonance in /i 
in combination with a plasmonic form, modified by the elec- 
tric cut-wire response of the SRR where necessary^**, which 
is essentially negative around the magnetic resonance fre- 
quency u!m- The effective impedance and index of refrac- 
tion of the HEM approximation for the first three unit cells 
of the LHM in propagation direction is shown in Fig. |9] as 
the colored points. In violation of the assumed effective 
medium picture we find again a cut-off of the negative reso- 
nant Re rihcm(w), this time at the lower edge of the first Bril- 
louin zone, Tiodgc = ~fccdgc/fc = — 7r/(fcL) where L is the 
length of the unit cell in propagation direction. The imaginary 
part of n does not have the expected form either. In Zhcm('jj) 
we expect two peaks for the LHM, one at Um and another 
one at the electric (effective) plasma frequency ujp which is 
the lowest frequency where e(cL') becomes positive. Though 
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the retrieved Zhcm(w) does show two such peaks, the position 
of the first one associated with ujm does not agree with the 
uj„i derived from nhem(w). This is the same issue as found 
for the SRR above. The corresponding effective permeabil- 
ity /ihcm('jj) and permittivity e\icm{^^) are shown in Fig. 1101 
The previously published resonance/anti-resonance coupling 
in the real parts around ti>.,„ is clearly visible, together with 
the appearance of the negative imaginary part in the permit- 
tivity and the misshapen absorption peak in the the perme- 
ability. Also for the LHM we can now confirm the length 
independence of the HEM approximation, apart from some 
additional noise, up to 10 unit cells in propagation direction. 
As for the SRR the PEM approximation of the LHM yields 
effective parameters free of all the artifacts seen in the HEM 
parameter which demonstrates again their origin in the peri- 
odic structure of the meta-material. In npcin(w) and Zpcm(ijj) 
we obtain an untruncated magnetic resonance and agreement 
of n and z upon the position of ljj„i- The permittivity Spcmi^) 
and permeability fipc„-,{uj) show a clean magnetic resonance 
with a symmetric absorption peak and the anticipated electric 
plasmonic behavior without any negative imaginary parts. In 
epcm(w) we can even recognize the beginning ascent of the 
imaginary part to the absorption peak at u; = contributed 
by the plasma resonance of the continuous wire. Note again 
the shift in uimp while the electric plasma frequency ui'^ is 
essentially the same for HEM and PEM approximation, al- 
though the descent to negative values for decreasing frequency 
is more rapid for epcm(ijj)- Also for the LHM these results 
are generic, ie. have been qualitatively reproduced for differ- 
ent sizes of the SRR and continuous wire components of the 
LHM and the unit cell. 



C. Continuous wire and cut-wire 

Not only the behavior of the SRR based meta-materials 
around the magnetic resonance but also that of the periodic 
arrangement of continuous wires or cut- wires are strongly af- 
fected by the periodicity. In the effective medium picture, the 
continuous wire is expected to expose a simple plasmonic per- 
mittivity going monotonously from negative to positive real 
part and crossing zero at a single plasma frequency cUp. In 
the HEM approximation of simulated continuous wire meta- 
materials^, eg. the isolated continuous wire array from the 
LHM discussed above, we observe the anticipated plasmonic 
behavior only at low frequencies up to the order of the plasma 
frequency cOp of the wire. At higher frequency we find a 
sequence of additional stopbands which can be explained as 
periodicity band gaps in the framework of the PEM model: 
whenever the residue class Re mod 27r/fc comes close 
to the edge of the first Brillouin zone a periodicity band gap is 
opened where a representative of Re nhcm('^) follows the line 
ficdgd^) and Im nhcm('jj) is significantly non-zero. At the 
boundaries of each of these band gaps the real and the imag- 
inary part of the effective impedance Zhom('jj) have either a 
zero or a pole which leads to the appearance of a alternating 
sequence of phony resonance-like structures in /ihoin(w) and 
ehcm(^)- The corresponding series of transmission above the 



"first" LUp have also been confirmed in experiments with thin 
metallic wires on PCB boards. Using the PEM approximation 
of the simulation data we can describe the scattering ampli- 
tudes in terms of a epcm(i-^) which does possess just the ex- 
pected plasmonic form in conjunction with an almost exactly 
/ipcm(^) = 1- The plasma frequency in ehem('i^) however 
does not coincide with the "lowest" ujp of the HEM approx- 
imation but appears moderately shifted to higher frequency. 
The HEM and PEM approximation for a meta-material com- 
prising a periodic array of continuous thin wires parallel to the 
electric field of the incident electromagnetic wave is shown in 
Figs. ^2 and [2] Both effective material approximations are 
virtually independent on the system length. 

The periodic array of cut-wires is of interest for two rea- 
sons: first, it can be used as a model of the electric response of 
the SRR and second, we could substitute the continuous wire 
in the LHM to tailor the collective £off (w) or to simplify its 
mechanical construction. For a sufficiently low resonance fre- 
quency uje the cut-wire shows a behavior analogue to the SRR 
but with the electric and magnetic parameters exchanged. We 
observe^"* all the previously described artifacts in the HEM 
approximation: the resonance/anti-resonance coupling, where 
the (electric) resonance appears this time in £hGm(^), accom- 
panied by an negative imaginary part in /iiicm(w), the cut-off 
of nhem(w) at the Brillouin zone edge, and so on. However, 
since in real meta-materials the cut- wire resonance usually ap- 
pears at much higher frequency, we observe an additional pe- 
riodicity band gap at lower frequency, well separated from the 
cut-wire response, qualitatively as shown in Fig.|3 The effec- 
tive parameters /ihomlt^) and ehcmlt^) at the lower boundary 
of this additional band gap will appear very similar to those 
of a cut- wire resonance at low frequency uje- Note that we 
can shift the cut-wire resonance to arbitrarily low frequency 
by reducing the gap in the longitudinal direction of the cut- 
wires, ie. making the finite wires almost as long as the width 
of the unit cell. The HEM and PEM approximation for such 
a meta-material of cut-wires parallel to the electric field of 
the incident electromagnetic wave is shown in Figs. 1131 and 
[T4l In either case, the effective parameters obtained from 
the PEM approximation do basically not show any artifacts 
around the cutwire resonance nor the periodicity band gaps 
and behave very much like expected from the naive effective 
medium picture. This gives us a criterion to distinguish the 
real cut-wire resonance in a meta-material from the possible 
phony one brought about by the periodicity. 



V. DISCUSSION 

When Veselago conceived the idea of heft-handed electro- 
magnetic material he essentially considered theoretical homo- 
geneous media for which there are no correspondents in na- 
ture. Though there are materials exposing a negative mag- 
netic response and other materials with a negative electric re- 
sponse, the challenge is to obtain both simultaneously in a 
material and, moreover, at an experimentally useful frequency 
in or above the microwave range. After Pendry proposed the 
first practical possibility to obtain controllable resonant mag- 
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FIG. 11: (color online) For the simulated continuous wire meta-material the effective index of refraction Weff (a;) and impedance ZcB{io) are 
shown. The colored curves represent the HEM approximation (Eqs. llOlll H of the simulation data for the first 3 unit cells, the solid black line 
the HEM(PEM) approximation, and the dashed line the PEM approximation (cf. Eqs. I51||521 obtained from the first unit cell data. 




FIG. 12: (color online) For the simulated continuous wire meta-material the effective permittivity Eoff (i^) and permeability fj,cs{i^) are shown. 
The colored curves represent the HEM approximation (cf. Egs. llOlfm of the simulation data for the first 3 unit cells, the solid black line the 
HEM(PEM) approximation, and the dashed line the PEM approximation (Eqs. I51||52> obtained from the first unit cell data. 
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netic and plasmonic negative electric response by the means of 
the geometric shape of electric conductors (SRR and contin- 
uous wire) in a periodic arrangement, numerical simulations 
have been conducted, attempting to verify the simultaneous 
negative magnetic and electric response. Note that because 
of the technical difficulty to obtain reliable complex reflec- 
tion amplitudes from experimental measurements numerical 
simulations are here the most important tool to address the 
question of negative electromagnetic response. The existence 
of negative refraction, which is merely a consequence of an 
anti-parallel phase and group velocity, has been demonstrated 
experimentally but provides no prove for the actual Veselago 
picture of simultaneously negative /i and e. For this purpose, 
at typical vacuum wavelengths of around ten times the size of 
the unit cell effective medium behavior of the meta-materials 
has been assumed throughout the literature such that an effec- 
tive permeability and permittivity could be obtained in a HEM 
approximation. Although the working frequency of a SRR- 
based LHM is theoretically arbitrary, the fabrication technique 
in experiments and the limited computer power in numerical 
simulations impose certain constraints on the size of the small- 
est structures (particularly the width of the gap in the SRR) 
in comparison to the size of the unit cell. As a consequence 
the ratio of the vacuum wavelength around the magnetic reso- 
nance to the size of the unit cell is confined as well. 

Continuum HEM approximation. The effective refractive 
index nhcm(t^), impedance Zhcm(^^) and, derived from those, 
permeability /Zhoin('^) and permittivity e\icm{'-^) published by 
usiiiiSaa and others2ii2iSLSa for the LHM, SRR and cut-wire 
meta-materials do only in first approximation meet the the an- 
ticipated effective response of the material. For the considered 
polarization (electric field parallel to the continuous wire and 
perpendicular to the gap bearing side of the SRR) and direc- 
tion of propagation in the SRR plane the resonant circulating 
currents inside the SRR ring should couple and respond only 
to the magnetic field, affecting the behavior of the effective 
The electric field couples to the continuous wire or cut- 
wire. However, it also couples to a separate electric resonance 
of the SRR caused by induced polarization currents oscillat- 
ing linearly in the gap-less sides of the SRR which are par- 
allel to the electric field. For other orientations, due to the 
anisotropy arising from the gap in the SRR ring, the electric 
field may also couple to the resonance of the circulating cur- 
rents making the total effective behavior a lot more compli- 
cated. In some cases violating the inversion symmetry of the 
unit cell also a magnetic coupling to the electric cut-wire re- 
sponse is possible. Though we obtain from the simulations 
a resonant magnetic (magnetic resonance frequency for 
SRR, LHM) and electric (cut-wire resonance frequency for 
SRR, cut-wire) response as well as basically the plasmonic re- 
sponse of the continuous wire and its corresponding contribu- 
tion in the LHM, the correspondence to the effective medium 
picture is spoiled by partially very significant anomalies: ii) 
Resonance/anti-resonance coupling. We expect the electric 
and magnetic response of the discussed meta-materials to be 
independent. However, whenever there is a resonance in Rc ji 
we simultaneously observe an anti-resonant behavior in Re e, 
and vice versa. The anti-resonant structures in the real part are 



accompanied by an negative imaginary part, (ii) Misshapen, 
truncated resonances. The divergence in Rc n appears to be 
cut-off at the edges of the first Brillouin zone and, in par- 
ticular, the negative regions of the magnetic resonance in ji 
and cut-wire resonance in e do not return from large nega- 
tive real part but seem to saturate in a rather shallow behav- 
ior. The corresponding absorption peak in the imaginary parts 
as misshapen and highly asymmetric too. {Hi) Discrepancy 
between n and z about the positions of the resonances. We 
expect the peaks (or zeros) in the index of refraction and the 
impedance to appear exactly at the resonance frequencies LOm 
and u)q, or the effective plasma frequency uj'p. From the simu- 
lations, however, we find different frequencies from n and 2, 
respectively. This lead for instance to an "internal structure" 
of the magnetic resonance as shown in Fig.|5] which can not 
easily be explained within the assumed effective medium pic- 
ture, {iv) Additional spectral structures. Apart from structures 
around the anticipated contributions of the meta-material's 
constituents ujm, and cj^ we observe a lot of additional 
structure, especially at higher frequency, which can not be ac- 
counted for The observed artifacts in the HEM approxima- 
tions are quite generic and have meanwhile been reported by 
different groups. An explanation for all these effects is given 
by the periodicity, see below. 

The occurrence of negative imaginary parts the permittiv- 
ity or permeability has been criticized by several authorsiii^. 
Indeed, a physical, passive homogeneous material may not 
possess any negative imaginary in /i and e. As long as the 
material is passive, it can only absorb energy from the elec- 
tromagnetic field. If there was any negative imaginary part 
in /i or £ one could device a geometry of this material which 
would violate the passivity. This requirement does not apply 
to the HEM approximation. The effective material is defined 
as the homogeneous material which reproduces the scattering 
amplitudes of the meta-material within the given geometry, 
ie. normal incidence to a homogeneous slab of finite length in 
propagation direction, if a length independent solution does 
exist. Our retrieved HEM parameters provide such a length 
independent description of the actual scattering amplitudes 
of the meta-material. For the (stationary) scattering of plain 
waves at a finite homogeneous slab in the continuum there are 
/i and e with negative imaginary parts which do not violate the 
passivity of the material, as long as much weaker conditions 
Im/i + Ime > A lm/i/|/i| + Ime/jel > are satisfied. 
If we further restrict the scattering setup requiring the thick- 
ness of slab to be an integral multiple of the unit cell, which is 
reasonable for this type of meta-material, we have even more 
freedom in jj. and e. The HEM approximation is only valid for 
a given geometry, there is not necessarily a physical material 
exposing the same material parameters in an arbitrary setup. 

Some authorsli have suggested a more general effective 
description^^ of the meta-materials which employ tensorial 
ji and e to take the coupling of the electric field to, and the 
anisotropy of the SRR into account. Though this is certainly 
an issue for arbitrary orientation of the SRR and will gain 
importance in more-dimensional materials, it is not directly 
related to the deviations from the effective medium behavior 
discussed above. 
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FIG. 13: (color online) For the simulated cut-wire meta-material the effective index of refraction ricff (cj) and impedance Zcft{uj) are shown. 
The colored curves represent the HEM approximation (Eqs. llOlllU of the simulation data for the first 3 unit cells, the solid black line the 
HEM(PEM) approximation, and the dashed line the PEM approximation (cf. Eqs. l51ll52l obtained from the first unit cell data. The weak 
additional structure in the PEM approximation close to the resonance is indicating the beginning breakdown of the approximation of the 
meta-material within our most simple periodic medium model. 
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FIG. 14: (color online) For the simulated cut-wire meta-material the effective permittivity Ecftii^) and permeability /icff (cj) are shown. The 
colored curves represent the HEM approximation (cf. Eos. 1101 fTTl of the simulation data for the first 3 unit cells, the solid black line the 
HEM(PEM) approximation, and the dashed line the PEM approximation fEQS. I51ll52> obtained from the first unit cell data. 
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Apart form some additional "noise", we could confirm the 
length independence of the HEM approximation for all con- 
sidered meta-materials up to 1 1 unit cells in propagation direc- 
tion. In contrast to previous work this became possible by the 
introduction of a explicitly isotropic material-discretization in 
the TMM and the use of a carefully chosen symmetry of the 
unit cell to avoid the occurrence of cross-polarization scatter- 
ing amplitudes'^. In the presence of cross-polarization terms 
the second polarization can contribute^" to the T and R of the 
considered polarization, eg. the transmission t22 ^ 1 for the 
passive polarization of the SRR which has the magnetic field 
in-plane and the electric field parallel to the gap-less sides can 
cut-off the transmission in the active polarization of the SRR 
in the stop-bands, tii{N) ~ + ^12 ^22^^ t2i + • • • , which 
should decay exponentially with the system length. This ef- 
fect can cause a phony length dependence in the HEM approx- 
imation which is diagonal in the polarizations. 

An other issue we have to address is the vacuum wavelength 
to unit cell size ratio. In the present paper we have: 1 (unit 
cell), 5 (SRR), 25 (low frequency SRR); a typical value for 
the meta-materials found in the literature is « 10 or worse. 
Strictly speaking we can expect effective medium behavior 
only in the cj ^ limit. As an approximation, it may hold 
if the wave length inside the structure is large compares with 
the typical length scale of the meta-material provided by the 
size of the unit cell. Two points are important here: The rel- 
evant wavelength is the wavelength inside the structure and 
can be much smaller than the corresponding vacuum wave- 
length if we approach the resonances at ijj,„ and uj^ and the 
magnitude of the real part of the effective index of refraction 
becomes large. Second, we cannot know a priori how large the 
above ratio has to be to reach a reasonable effective medium 
behavior Our simulations support both points as the devia- 
tions from the expected behavior happen essentially when the 
wavevector inside the structure, q = n{u)) k, becomes compa- 
rable with the edge of the Brillouin zone, ie. the wavelength 
comparable to the unit cell size; for the low frequency SRR 
we find pretty good effective medium behavior without the 
afore mentioned anomalies apart from a small region around 
the resonance. The simulations indicate that to obtain a rea- 
sonable effective medium behavior the wavelength to unit cell 
size ratio has to be of the order of 30. The less losses occur in 
the resonances the larger the peak in n, and consequently the 
required wavelength to unit cell size ratio will be. 

Continuum PEM model. The analytic HEM approxima- 
tion of the periodic medium model proves that the observed 
deviations from the anticipated effective medium behavior in 
the HEM approximation of the real simulated meta-materials 
can be caused by the periodicity or, more precisely, the re- 
duced translational symmetry to the discrete group generated 
by the unit cell. In particular, the model shows that those ar- 
tifacts are not related to the actual geometric definition of the 
SRR and LHM resonances since the homogeneous core in the 
model just uses the resonant forms of Ijl{lo) and e((jj) without 
any reference to their microscopic source. 

The artifacts observed in the HEM approximation ulti- 
mately originate from the occurrence of band gaps introduced 
by the periodicity, analogously to, for instance, the band- 



structure in a crystal. In our case, the only substantial differ- 
ence is the explicit frequency dependence of the microscopic 
material properties, ie. the resonant forms of /i(ijj) and e(aj) 
for the homogeneous core of the periodic medium model's 
unit cell. These periodicity band gaps are distinct from the in- 
trinsic band gaps which arise directly from the negative prod- 
uct e/i in the emulated response of the constituents. To illus- 
trate the behavior in the periodicity band gaps we assume for 
the moment ^{uo) and e{Lo) of the core to be real. Then n{Ld) 
and z{uj) in Eq. I^^are simultaneously either real or imag- 
inary and consequently the r.h.s. of ( I23> is real as well. A 
periodicity band gap occurs whenever the rh.s. grows out- 
side the domain [—1,1] of the cosine for a real argument: 
For the principal branch^^ UeS acquires a non-zero imagi- 
nary part inside the gap and Re nofi (^) retains the value of 
zero or tt/ [kL) for all rh.s. of equation (I23> above 1 or be- 
low — 1, respectively. Adding a small imaginary part m ^{lu) 
and e{Lo) of the core, as should be a good approximation for 
the emulated meta-material if we do not come too close to the 
resonances (absorption peaks), adds a small imaginary part 
to ricff(w) and causes Re ncff(aj) to deviate slightly from 
or tt/ {kL) towards tt/ (2fcL) inside the periodicity band gap. 
Having thus established the confinement of the effective in- 
dex of refraction. Re noff(w), to the edge(s) of the fist Bril- 
louin zone or to zero, the coupling of /i and e follows as a 
direct consequence^'. In simple words: If = efi is con- 
fined and either one of e, /i exposes resonant behavior the 
other has to go to zero simultaneously. This also explains why 
the resonance/anti-resonance coupling and the negative imag- 
inary parts occur only across the periodicity band gaps asso- 
ciated with the resonances but not outside. For the general 
case this qualitative behavior is complicated by the non-zero 
imaginary parts. 

Note again that we can obtain the real part of the index of 
refraction only as a residue class Re Tice{u!) mod 27r/(fcL) 
which becomes immediately clear either from the length- 
independence of the HEM approximation of the PEM accord- 
ing to equations ( IT9l I20> or from the argument about the 
simultaneous congruences for different system lengths dis- 
cussed in sectionin] Therefore the pieces of Re nhcm('jj) that 
follow the multiples of the first Brillouin zone's edge just co- 
incide with either zero or the first Brillouin zone's edge (upper 
and lower are equivalent) itself. 

An other important observation is that a periodicity band 
gap may occur in between the magnetic resonance uj„i and 
the cut- wire resonance uj^- Unfortunately, the corresponding 
effective e looks like an electric resonance with the attendant 
anti-resonant structure in /i. This may easily be mistaken as 
the cut-wire resonance which would be expected to follow the 
magnetic resonance as the next feature in the frequency spec- 
trum. Only in the low-frequency limit the latter behavior is 
actually observed. Not that the parameter dependence of the 
phony electric resonance frequency at the lower edge of the 
periodicity band gap will qualitatively resemble the behavior 
of the real cut- wire resonance frequency Ue. 

In comparison to the TMM simulations we see, however, 
too much structure, ie. a series of periodicity band gaps in- 
stead of only one, around the intrinsic resonances. The basic 
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difference between the TMM simulation and the analytic cal- 
culation of the scattering amplitudes in the continuum is the 
presence of a finite discretization mesh in the TMM, which 
implies a smallest distance, hence in turn a largest supported 
momentum. Obviously, this hmitation will become visible 
where ricei^^) grows large. 

Lattice PEM approximation. In our TMM simulations the 
lattice version of the PEM does better correspond to the nu- 
merical data than the continuum PEM. This does also apply to 
independent Microwave Studio simulations which, in contrast 
to our TMM simulations, utilize a non-uniform discretization 
of the meta-material. If the discretization mesh is chosen finer 
the lattice PEM gradually approached the continuum PEM be- 
havior. 

As discussed above the HEM and the PEM approximation 
of the real meta-material are basically length-independent, 
longer systems expose the same spectral features as the first 
unit cell. However, systems that consist of more than a sin- 
gle unit cell in propagation direction do contribute additional 
tiny resonance-like structures in the effective material con- 
stants /ihcm(^^), Ehcmlw) which are also present in the ef- 
fective parameters of the PEM approximation, /ipcm(^) and 
£poin(w) (not shown). For a slab of N unit cells in prop- 
agation direction these additional structures appear as tiny 
band gaps, quite similar to the periodicity band gap discussed 
above, at frequencies where 7T,pem(^) ~ rrm / {NkL) with 
m e [~N,N] C Z, ie. whenever the effective refractive index 
derived from the simple unit cell comes close to the multiples 
of the first Brillouin zone's edge for the whole slab. The ad- 
ditional structures weaken and eventually ceases to be visible 
in the low-frequency limit as can be seen, for instance, for the 
SRR in Figs. [7] and |S] Note that again the behavior is generic: 
it appears in our TMM simulations for LHM, SRR and also 
cut-wire and continuous wire meta-materials (not shown), it 
has also been verified in Microwave Studio simulations using 
a different numeric technique. As there is no such length de- 
pendence in the analytic periodic medium model, we interpret 
this additional "noise" as a limitation of the PEM approxima- 
tion of the real meta-material which starts to see some inter- 
nal structure apart from the explicit periodicity. It is not yet 
clear whether the groups of N peaks directly at the magnetic 
resonance are related to this problem. For these there are at 
least two other interpretations: They may be caused by the 
coupling of successive SRRs in propagation direction which 
would lead to the splitting of the resonance frequency a;,„ as 
for the eigenfrequency of coupled identical oscillators, or the 
finite accuracy of the numeric simulation data could, in partic- 
ular around the resonances, lead to an residual explicit length 
dependence. 

We demonstrate above that the PEM approximation of real 
LHM and SRR meta-materials is good in the region around 
bJm if the corresponding vacuum wavelength in as small as 
five times the size of the unit cell. Then we obtained the an- 
ticipated effective behavior in /ihcm('^), £hcm('^) and all the 
anomalies and additional features seen in the HEM approxi- 
mation arose from the explicit periodicity. If we move to even 
higher frequency where the vacuum wavelength is close to 
the size of the unit cell, also /ihom(w) and ehcm(w) start to de- 



velop unexpected features like additional magnetic response 
around the cut-wire resonance. Although even then the PEM 
spectrum is not nearly as erratic as the corresponding HEM 
spectrum, this indicates the breakdown of the PEM approxi- 
mation. This is not surprising since we now reach the photonic 
crystal limit and the internal structure of the meta-material' s 
constituents must become visible and no "effective descrip- 
tion" should be possible anymore. 

Momentum dependent parameters. In order to take into 
account the periodic structure of real meta-materials, we con- 
sidered the PEM model as a very simple explicit example of 
a periodic medium. Instead of qualifying the specific geom- 
etry of the model used, we may alternatively introduce k- 
dependent effective parameters fj,{k,LLj), e{k,uj) to describe 
the spatial distribution of the electromagnetic response in 
the PEM. These /c-dependent effective parameters completely 
characterize the effective medium model. For the PEM model 
defined in Fig.[l]we find by Fourier transformation e{k, oj) — 
(27r)^^/^ / dz e{z, uj) e''^^ the representation 

eik,u;) = (53) 



27r ecorc(t^) 



ikL 



J2s{k 



27rTO 



and correspondingly for ij,{k,uj). As a generalization of the 
periodic medium approximation we could further ask, which 
arbitrary /s-dependence of e{k, uj) and fi{k, uj), ie. which spa- 
tial distribution of the effective material parameters, describes 
a given meta-material best. Although this might be desirable 
as a descriptive tool for the engineering of meta-materials, it 
is clearly beyond the scope of this paper. 

What is the actual left-handed band? Finally we want to 
comment on the actual extent of the left-handed interval of the 
LHM as it concerns experiments and applications. Obviously, 
the bands with Re ries < retrieved via HEM and PEM ap- 
proximation differ considerably in width (Fig.|9jl which raises 
the question, where to expect the left-handed behavior. We 
argue, that the correct region is given by the HEM approxi- 
mation. If there is a length-independent HEM approximation, 
the scattering behavior of the meta-material can be described 
assuming plain-wave solutions inside the homogeneous unit 
cell. This plainwaves will possess a wavevector q related 
to the vacuum-wavevector k by the retrieved negative index 
of refraction, q = nhcm(w) k. They coincide with the non- 
periodic factor of the Bloch-waves describing the periodic 
medium, ie. coincide with the Bloch-waves at the edges of the 
unit cells in the meta-material. Therefore a (usually damped) 
plainwave with negative phase-velocity will exist inside the 
meta-material whenever Re rihom < 0. This interpretation is 
also supported by experimental measurements'^. Note that this 
frequency interval is wider than the interval with simultane- 
ously negative Re ehcm and Re /ihcm- The PEM approxima- 
tion indicates that the isolated local response of the SRR and 
wire, without the effects of periodicity, would lead to a much 
smaller left-handed band (compare the behavior of the low- 
frequency SRR, Figs.0|SJl. The modifications of the generic 
response of SRR and wire by the band-structure, in partic- 
ular by the emergence of periodicity band-gaps, arising from 
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the inherent periodicity of the meta-material, greatly enhances 
the width of the negative index band in meta-materials which 
see strong artifacts from the periodicity. This concerns virtu- 
ally all published meta-materials with a vacuum-wavelength 
to unit cell length ratio smaller than approximately ten. At 
much lower frequency, truly effective homogeneous behavior 
will emerge, the periodicity band-gaps disappear, and HEM 
and PEM description coincide. 

Geometry of the PEM model. As we demonstrated in 
section |ni| for the quasi one-dimensional scattering prob- 
lem of a system comprised of an integral number of unit 
cells in propagation direction, we can exactly describe any 
HEM by a family of PEM, parametrized by the geometry 
Q = (ria, nd^ rib) of the periodic medium model (cf. Fig.[0, 
in terms of effective parameters /ipem[e](a;) and £pcm[e]('^)- 
For a simulated meta-material, in general all effective pa- 
rameters, /iiiom(w) and £hcm(w) as well as /ipoin[g;] (w) and 
£pcin[g] (t^)' will show "unphysical" behavior caused by the 
internal spatial structure of the unit cell and the periodic- 
ity of the meta-materials. If for a real meta-material the lo- 
cal electromagnetic behavior of the constituents can be ab- 
stracted from their geometrical form, approximated by sim- 
ple resonant and plasmonic response functions /isRR(ijj) and 
Ev^irciuj), respectively, and separated from the effects of the 
periodicity (ie. band-structure), then there is a particular ge- 
ometry of the PEM which approximates physical parameters 
A*pom[e](t^) ~ A^srrC^^) and epe,„[g](w) w ewiro(w) without 
the usual artifacts discussed in this paper In our simulations 
the best such approximation was obtained if the lattice PEM 
containing a single plain of scatterers in the unit cell, ie. for a 
geometry G — (5,1,4). 

VI. CONCLUSION 

We have investigated the influence of the inherent peri- 
odic structure always present in meta-materials which are 
build from the repetition of a single unit cell on the effec- 
tive medium approximation. It has been shown analytically 
that all the previously observed violations of the anticipated 
effective medium behavior of the (single-ring) SRR and LHM 
involving a single magnetic and a single electric resonance 
can be explained in term of the periodic structure: A very 
simple stratified periodic medium model involving slabs of 
vacuum alternating with slabs of a homogeneous material 
with simple resonant and e(cj) can reproduce all the 

artifacts like resonance/anti-resonance coupling in /ihcm('-^) 
and ehem(w), negative imaginary parts in either e or /i, trun- 
cated, misshapen resonances, additional band gaps and the 
complicated high-frequency behavior found in the HEM ap- 
proximation of numerically simulated SRR arrays and LHM, 
but also meta-materials build of continuous wires and of cut- 
wires. In good approximation, the effective behavior can be 



decomposed into a effective behavior of the constituents of 
the meta-material and an explicit contribution of the periodic- 
ity. Remarkably, the average contribution of the constituents 
like the single split-ring behaves much like expected from the 
assumed homogeneous medium picture, which can only be 
justified in the low-frequency limit where the wavelength in- 
side the structure is large compared to its geometrical size, up 
to frequencies where the vacuum wave length becomes com- 
parable to the size of the unit cell. This allows a more reliable 
effective description and interpretation of real meta-materials 
in terms of a periodic effective medium (PEM) instead of the 
conventional homogeneous effective medium (HEM) with all 
the hard to understand artifacts. The effects caused by the pe- 
riodicity are generic, they do qualitatively not depend on the 
particular geometry chosen for the meta-material and univer- 
sally apply to SRR, LHM, continuous wire and cut-wire ma- 
terials. Obviously, the impact of the periodicity is noticeable 
throughout the range of the Xq/L ratio « 5 ... 10 (ie. vac- 
uum wave length / unit cell length) actually found in published 
simulations and experiments for left-handed and related meta- 
materials. Our simulations indicate that an unencumbered ho- 
mogeneous effective medium behavior, though reachable in 
the low-frequency limit, would require a \q/L ratio in the or- 
der of 30 or larger which is geometrically not easy to obtain in 
real samples. We investigated the difference between the con- 
tinuum and the lattice formulation of the PEM approximation 
and found the latter to be better suited for application to our 
numerically simulated scattering data for real meta-materials 
obtained with a lattice-TMM implementation. 

The PEM approximation may provide an valuable tool to 
understand the various features observed in the scattering 
spectra of real meta-materials in experiments and simulations. 
In the present paper we essentially discussed the vicinity of 
the magnetic resonance of the SRR. Further work shall em- 
phasize on the frequencies above the magnetic resonance Um 
including the cut-wire resonance uj^- 

We expect the impact of the meta-material's periodicity to 
be noticeable also in higher dimensional structures. Because 
the unit cells of those structures tend to be more complicated 
including couplings of SRRs in the different directions the 
separation of the "real" effective response of the constituents 
and the structures produced by the periodicity constitutes an 
even more immanent issue for understanding. 
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We consider the transfermatrix T for a planar stratified medium in 
z-direction. Because div D = 0, div B = 0, the source-free 
Maxwell equations only allow for four independent field com- 
ponents: F = {E^,Ey,H^,Hy), ¥{z + dz) = T{z)W{z). 
The eigensystem {k,¥k} of the transfermatrix inside the strata. 



l{z)Wk{z) 



''¥k{z), defines the plainwave base (includ- 



ing propagating and evanescent waves) for the scattering problem. 
Local 2;-isotropy implies that with each k also — fc is eigenvalue 
such that the wave basis splits into right- and left-going waves, 
while translational invariance of the stratified medium ±2 fac- 
torizes the transfermatrix with respect to the momentum fcy par- 
allel to the strata which is preserved across interfaces. Within a 
fc|| )-subspace the E and H fields are dependent, coupled by 
the Maxwell equations, such that we can decompose F as a lin- 
ear combination of only two polarizations, eg. TE and TM mode. 
For the TE mode, the electric field E is parallel to the interface 
and we can express it inside each stratum as a linear combination 



+ Be~ 



. For the TM mode the same applies to 



the magnetic field H. In either case we can write a transfermatrix 
for the wave amplitudes across a single interface as {A' B')^ — 
%y{z){A Bf, where Tkk'{z) = </)-/(z) T^^, (/.^(z) can be 
written as product over phase factors (f) and the location indepen- 
dent matrix Tfej./ which coincides with the transfermatrix for the 
fields (V'fc V'-fc) = (^e'''^ Be"*''^) instead of the ampfitudes. 
In explicit form we have 
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with Ckk' = p.'k/{fik') for the TE and e'k/{ek') for the TM 
mode. The transfermatrix for a vacuum- terminated stack of lay- 
ers has the composition-property, T13 — T23T12, and inserting a 



virtual vacuum interface Tk^k„^i = T^^^^k^^ 
each adjacent two layers we get an expression 



Ttotal {k 
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W rslab(z„ - Zn-l) 
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0vac(2o) 



with the interfaces located at Zn and the local transfermatrix of a 
single layer rsiab(d) = Tf,^k^^„ ^k„{d) Tk^^^k„ which is given 
explicitly in SectionllTI 
"'^ We define the principal branch of a,rg{z) as —n < arg z < n 
with a branch-cut along the negative real axis, continuous from 
the second quadrant. Then the principal branches of ^/z = 
\z\^^^ exp(i arg(2) /2) and log (2) = log \ z\+i arg(2) can be de- 
fined with the same branch-cut and continuity from above. Using 
2cos(2) — exp(j2) + exp(— 12) we obtain the principal branch 
arccos(2) = —i log{z + i\J\ — z^) with a branch-cut on the real 
axis for 2 G (1, 00), continuous from below, and 2 G (—00, — 1), 
continuous from above. 



